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Abstract 

^ . We investigate the influence of the neutron halo and the breakup channel 

in ^He + ^ss-q f^gjon at near-barrier energies. To include static effects of 
the 2n-halo in ^He nuclei, we use a single-folding potential obtained from an 
^ ■ appropriate nucleon-^^^U interaction and a realistic ^He density. Dynamical 

■ effects arising from the breakup process are then included through coupled- 

^ . channel calculations. These calculations suggest that static effects dominate 

O ! the cross section at energies above the Coulomb barrier, while the sub-barrier 

2 I fusion cross section appears to be determined by coupling to the breakup 

channel. This last conclusion is uncertain due to the procedure employed to 
measure the fusion cross-section. 
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1 Introduction 



The recent availability of radioactive beams has made possible to study re- 
I actions involving unstable nuclei [1] . Several of the light neutron and proton 

rich nuclei exhibit halo structures, with a compact core plus one or two loosely 
bound nucleons. For example, ^^Li and ^He are two-neutron, borromean halo 
nuclei, while ^^Be and ^^C are one-neutron halo nuclei. The isotope *B has 
been confirmed to have a one-proton halo, while ^^F is a normal nucleus in 
its ground state but becomes a one-proton halo in its first excited state. 

Reactions induced by these nuclei are important in processes of astrophys- 
ical interest, among others. We ask the question of how the above systems 
fuse, in particular how the fusion induced by these nuclear species behaves 
as a function of bombarding energy, especially near the Coulomb barrier. 
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The main new ingredient in reactions induced by unstable projectiles is 
the strong influence of the breakup channel. One important feature of these 
loosely bound systems is that they exhibit the so-called soft giant resonances 
(pygmy resonances), the most notorious of which is the soft dipole resonance, 
very nicely confirmed in ®He by Nakayama et al. [2]. 

In the case of not too unstable projectiles, the effect of the breakup chan- 
nel in the fusion cross section at sub-barrier energies is, as in the case of 
stable beams, to enhance it. At slightly above-barrier energies, however, the 
situation is qualitatively different from the case where only stable nuclei are 
involved. The contribution from the breakup channel to the fusion reaction 
is strongly inffuenced by the low probability that all fragments are captured. 
Thus, in this case, the fusion cross section is partitioned into a complete and 
one or more incomplete fusion contributions. 

Recently, nuclear reactions involving the neutron-rich nucleus ^He have 
attracted considerable attention [3]. In particular, very interesting exper- 
imental data on the fusion of He isotopes with ^ss-y j-^g^yg been obtained 
[4]. These data show an enhancement of several orders of magnitude of the 
^He-|-^^^U fusion cross section with respect to that of ^He-|-^^^U. The phys- 
ical process leading to this result has not yet been established. A natural 
candidate is the coupling with the breakup channel. This led us to develop 
a simple theoretical model to estimate statical and dynamical effects of the 
breakup channel on the complete and incomplete fusion cross section in the 
^He-|-^^^U collision. The extension of the model to study fusion induced by 
other radioactive beams is straightforward. 

The paper is organized as follows. Section 2 describes the calculation 
of the static effects brought about by the presence of a nuclear halo. The 
coupling to the breakup channel is performed, by means of schematic coupled- 
channels calculations, in section 3. Our conclusions are presented in the last 
section. 

2 Static effects from the 2n-halo 

The weakly bound neutrons in ^He are expected to inffucncc the fusion cross 
section in two ways. Firstly, by the static effect of barrier lowering due to the 
existence of a halo. Secondly through the coupling with the breakup channel. 
In this section we consider the first of these effects. 
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Owing to the two weakly bound neutrons in ^He, the nuclear density has 
a long-range tail and so does the real part of the optical potential describing 
the ^He-target collision. In this way, the potential barrier is lowered and 
the fusion cross section is enhanced. In order to account for this effect, we 
use a single folding model do describe the real part of the nuclear ^He-^^^U 
interaction. This potential is given by the expression 

14r(r)= J VnMr - r') p{r') dh' . (1) 

Above, f „_r(r — r') is a phenomenological interaction between a nucleon and 
the ^^^U target nucleus and p(r') is a realistic ^He density, containing the 
contribution from the halo. The nucleon-^^^U interaction is obtained from 
studies of the collision of low-energy neutrons with heavy target nuclei in the 
actinide region. It can be written (discarding the spin-orbit part) [5] 

Vn-T{x) = -Vofr{x), (2) 

with 



AT" y 

Vo = 1 50.378 - 27.073 ( ^ j - 0.354 ElJ (Mel/) (3) 
and 

" l + exp[{x-Rr)/ar] ' 

with the parameters Rr = 1.264 A^''^ fm and = 0.612 fm. The total optical 
potential is then given by 

U{r) = VN{r) + Vc{r) - iW{r) . (5) 
Above, Vc(r) is the usual Coulomb interaction in nuclear collisions. 

Voir) = r > Rc ^ 1.2 {A]/' + A'/') 

r < Rc, (6) 
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and W{r) is a volumetric strong absorption potential with small values for 
both its radius and diffusivity. We use the parametrization 

W{r)^Wo Mr), (7) 

with Wo — 50 MeV and /j(r) a Wood-Saxon shape as in eq.(4) with 

Ri = 1.0 (Ap/^ + A^^^) fm, ai = 0.10 fm. (8) 

As a test, we apphed the above procedure to '^He + ^^^U fusion. The 
nuclear potential was evaluated by eq.(l) using a Gaussian density. We write 

p{r) = C exp (-r^'j^) (9) 

and choose the parameters C and 7 as to give the correct normalization and 
experimental r.m.s. radius. That is 

Jpir)d'r^A; J p{r)d'r ^ r'^^, . (10) 

In the present case, we set A = 4 and Vrms = 1-49 fm [6]. The fusion cross 
section obtained with our optical model calculation with the single folding 
potential is shown in figure 1 (thin solid line), in comparison with the data 
of Trotta et al. [4] and the data of Viola and Sikkland [7] . The agreement is 
very good. Since the calculation contains no free parameter, this agreement 
indicates that the procedure is reasonable. 

We now consider ^He + ^^^U fusion. Firstly, we disregard the existence 
of the ^He halo and repeat the above procedure. We parametrize the density 
as in eq.(9) and scale the density and r.m.s. radius to ^He. That is, we set 
in eq.(lO) A — 6 and Vrms — (6/4)^^^ x 1.49 fm. This density is then used 
in eq.(l) and the folding potential is determined. The fusion cross section 
calculated with this potential is shown in figure 2 (dashed line), in compar- 
ison with the data [4]. The agreement is poor throughout the considered 
energy range. We now take into account the existence of the ^He halo, re- 
placing the gaussian of eq.(9) by a reahstic parametrization [6] of the ^He 
density, based on the symmetrized Fermi distribution of ref. [8]. It leads to 
the r.m.s. radius Vrms = 2.30 fm. Using this density in eq.(l), we obtain a 
potential which includes contributions from the ^He-core and also from the 
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2n-halo. The resulting fusion cross section is represented by a solid line in 
figure 2. We note that the agreement with the data at above barrier ener- 
gies {Ecm. > Vb ~ 22.3 MeV) is considerably improved. Since the Coulomb 
barrier height is reduced by the attractive contribution from the halo, the 
cross section becomes larger. However, at sub-barrier energies the agreement 
remains very poor. The theoretical prediction for the fusion cross section is 
still several orders of magnitude smaller than the experimental data. 

3 Coupled channel effects 

It is well known that the coupling between channels enhances the fusion 
cross section at sub-barrier energies [9]. Therefore, coupled-channel effects 
should be taken into account in a theoretical description of the fusion process. 
However, in the case of coupling to the breakup channel the situation is 
more complicated since the breakup channel involves an infinite number of 
continuum states. A possible treatment of the problem, used in refs. [10, 11], 
is to use continuum discretization to reduce it to a finite number of channels. 
The situation is still more complicated in the breakup of ^He, since it breaks 
up into three particles. In the present work we schematically replace the 
breakup channel by an effective channel [12]. This state has energy equal to 
the breakup threshold and carries the full strength of the continuum. This 
procedure is justified when breakup occurs through a low-lying long-lived 
resonance (with a half life much larger than the collision time), as it seems 
to be the case with ^He [2]. Since the kinetic energy of the relative motion 
between the ^He-core and the neutron pair is neglected, this approximation 
tends to overestimate the importance of the coupling to the breakup channel. 
Therefore the simplified model of the present work should provide an upper 
limit for the fusion cross section. 

The starting point of the coupled channel method is the Schrodinger 
equation for the colliding system, 

i/*(r,0=£^*(r,0, (11) 

where r is the projectile-target vector, ^ stand for the relevant intrinsic co- 
ordinates, E is the total energy in the center de mass frame and H is the 
total Hamiltonian of the system . One then performs the channel expansion 
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of the wave function 

*(r,O = EV'a(r)0a(O, (12) 

a 

where (t>a{C) denotes an intrinsic state with energy and '0a (r) is the relative 
motion wave function in channel-a. Substituting this expansion in eq.(ll), 
we obtain the coupled- channel equations 

{E^ - H^) M^) = E^a/3(r) M^)- (13) 

Above, Ea = E — Sa and Ha — T+Ua{r), where Ua{r) is the optical potential 
in channel-a. The channels are coupled through an interaction V(r,^), with 
matrix-elements in channel space given by 

uap{r)^ |de€(0^(r, 0/3(0- (14) 

For practical purposes, it is convenient to carry out angular momentum ex- 
pansions. The wave function is then written as (see e.g. [13]) 

*W(aoz^oko; r) = i^M MM - u,) huo) >^MM-.o)(ko) 

X Eyai^io (15) 

al ^0?^ 

and using this expansion in eq.(ll) one obtains the angular momentum pro- 
jected coupled channel equations 

^J2^al,a'l'ir) Ua'l',Oloika',r) . 
a'l' 

In the present calculation, a takes only the values (elastic channel) and 
1 (effective breakup channel). For the energy of the breakup channel we 
used ei = 0.975 MeV, which corresponds to the breakup energy. As said 
above, this means we neglect the kinetic energy of the relative motion of the 
fragments after breakup. 
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We initially consider the coupling interaction as the electric dipole term 
in the multipole expansion of the electromagnetic interaction between the 
projectile and the target. This is based on the idea that, in order to break 
a very weakly bound nucleus, only a small perturbation is needed. The fact 
that the breakup cross section for those nuclei is very large, suggests that 
this process is important. 

In the case of a electric dipole interaction, the coupling matrix elements 
are [13] 




with 



(17) 



A^eZT^B{El,Q^l) {-Y+^ (18) 



Above, ^ Q Q j and | ^ g | the usual 3J and 6J symbols 

[14]. Note that the above matrix-elements are fully determined, except for 
the value of the reduced transition probability B{E1, 0-^1). 

Solving the coupled channel equations, one obtains the fusion cross section 
by the formula^ 



= {2nf ^ E (^i^^ W» ^'^') ■ (19) 

^ a=0 

The method of the present work was used to evaluate the fusion cross 
section in the ^He+^^^U coUision. We used the optical potential discussed in 
the previous section, which include the static effects of the halo. The coupling 
matrix-elements were given by eqs.(17) and (18), with the B[E1, 0^1) given 
by the cluster model [1], 

3;?2p2 

BiEl,0^1) = —^ • (20) 

lb7rei//2n-4He 



^The constant (27r)^ in the expression for ap arises from the normalization factor 
(27r)"^/^ adopted for vi'^^ 
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Above, ei is the energy binding the dineutron to "^He in the ^He nucleus and 
/i.2n-4He is the Corresponding reduced mass. Taking the numerical value of 
eq.(20), we obtain B{El, 0^1) = 1.37 fm^. 

Recently Hagino et al. [10] have shown that the effects of the nuclear 
coupling may extend quite far in the case of weakly bound nuclei. In order to 
estimate the additional dynamic effects arising from the nuclear interaction, 
we must include the coupling due to the nuclear potential. Since we use an 
effective channel to describe breakup states, the calculation of the nuclear 
form factor is a complicated task. For the estimates of the present work, we 
considered the nuclear interaction potential associated to ®He breakup to be 
the difference between the sum of the nuclear potentials between ^ss-jj ^^^^ 
^He and the dineutron, and between ^^^U and the ^He projectile, i.e. 

y,^,(r, x) = l^4He(r + x/3) + ^„(r-2x/3) - l^6He(r) . (21) 

Above, X is the vector going from the di-neutron to the ^He cluster, V2n is 
twice the potential of eq.(2) and V4hc, and Ve-^e ^re the folding potentials of 
the previous section. We carry out the angular momentum expansion 

Ct = E^A.(r)n;(x) V^M (22) 

and keep only the dipole term (A = 1). In this way, the nuclear form factor 
is 

Fj^^,(k;r) = / drM^) V,''{r,x) u,(k,x) , (23) 
Jo 

where (f)o{x) is the radial function associated to the bound state of the 2n-^He 
system and Ui{k,x) the / = 1 continuum radial wavefunction for the same 
system, with energy Ef^ = /i^/c^/2/X2n-4He- Both functions are calculated 
using the radial Schrodinger equation associated to the internal coordinate 
X. The depth of the V^n-^He potential was set in order to have the second 
S-state with energy Eq = —0.975 McV (to be consistent with Pauli Principle 
we discarded the first S-state). Owing to the normalization of ui{k,x), the 
above form factor vanishes in the A; — > oo limit. However, the absolute 
strength of Fj^^^ should be treated as a free parameter, since the final state 
is an effective channel. In this way, we adopt the form factor 

(r) = F, /(r) , (24) 
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with 



/(r) = lim 

fe— >oo 




(25) 



To estimate the strength Fq, we adopt the following procedure. Firstly, we 
evaluate the Coulomb form factor as we evaluated the nuclear one. Instead 
of using B{E1,0 — >• 1) = 0.59 fm^, we calculate reduced matrix elements 
of the dipole term in the Coulomb coupling using the analog of eq.(23). The 
resulting Coulomb and nuclear dipole form factors are shown in figure 3. 
Since the dipole term of the nuclear coupling cannot be written as a product 
of a function of r times a function of x, as can the Coulomb coupling to a good 
approximation, the shape of the nuclear form factor depends on the energy of 
the continuum state in the x-space. However, the shape of the nuclear form 
factor does not change much as — > 0. Although both form factors go to 
zero in this limit, they decrease by a common factor. In figure 3, we show the 
Coulomb and the nuclear form factors for a very low energy in the continuum. 
We see that the ratio of these form factors changes strongly with the radial 
distance. The Coulomb form factor dominates at large separations while the 
nuclear form factor is larger at small separations. They have approximately 
the same strength at r ~ 16 fm. In the present calculation, we use the 
experimental B[E1, 0^1) value and choose the parameter Fq such that the 
ratio between the nuclear and the Coulomb form factors is maintained. 

Figure 4 shows the ^He + ^^^U total fusion data in comparison to the static 
(dashed line) calculation of the previous section, and two coupled channels 
calculations. The thin line is the coupled channel calculation restricted to 
Coulomb breakup. We notice that the cross section at high energies is little 
affected by the inclusion of the breakup channel. Although the sub-barrier 
cross section is larger than that found in the previous section, it remains 
much smaller than the experimental values. 

The solid line is the calculation including also the nuclear coupling. We 
notice that it also changes little the cross section at high energies, and al- 
though the nuclear coupling affects more the fusion cross section at sub- 
barrier energies, the slope remains much larger than that suggested by the 
data. Changing the strength or diffuseness parameters of this couphng does 
not change this behavior. 

It should be pointed out that the coupling with excited states of ^38^ jg 
not likely to be relevant for this issue, since they were not necessary for the 
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description of the ^He+^^^U fusion data, considered in section 2. As our cal- 
culation should provide an upper limit for the cross section, the experimental 
fusion cross section at the lowest energies cannot be explained through our 
calculations. However, one should keep in mind that in the calculations pre- 
sented here we have not included effects due to coupling to other channels 
other than breakup, and in particular the transfer channels. As transfer 
close to the optimal Q-value may be quite important at sub-barrier energies 
[15], coupling to those channels, which should not affect much the ^He-|-^^^U 
fusion, is expected to influence strongly sub-barrier ^He-|-^^^U fusion. This 
could also be the case for the ^He+^''^Bi total fusion cross section, where the 
data [16] show a similar trend as sub-barrier energies. 

4 Conclusions 

We have investigated static and dynamic effects on the ^He-|-^^^U fusion cross 
section. Static effects of the halo were taken into account through the use of 
an appropriate optical potential. This potential was obtained by the single 
folding model, with a nucleon-targct interaction which is able to reproduce 
the ^He-|-^^^U fusion cross section and from a realistic ^He density. Dynamical 
effects were considered in a simplified coupled channel calculation, in which 
the breakup channel was represented by a single state with energy ei = 0.975 
MeV (the threshold for ^He breakup), concentrating all the low energy dipole 
strength. From our calculations we concluded that the static effects dominate 
the behavior of the fusion cross section at energies above the Coulomb barrier. 
The dynamic coupling to the breakup channel is important mostly below the 
barrier. It may be separated into the Coulomb and nuclear contributions. 
Although the breakup process takes place at large distances, we have shown 
that the coupling with the breakup channel cannot reproduce the main trends 
of the data in the sub-barrier region. We believe that demonstrates that a 
full description of the ^He-|-^^^U fusion cross section at sub-barrier energies 
requires the inclusion of neutron-transfer channels. We point out that a 
similar enhancement of the sub-barrier fusion cross section has also been 
observed in the collision of ^Hc with ^'^^Bi. 

After the completion of this paper we have learned [17] that the data of 
Trotta et al. [4] have been reanalized and new data with a different experi- 
mental set up have been taken. The new set of data seems to indicate that 
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the large enhancement at sub-barrier energies is due to transfer-fission, rather 
than fusion-fission events. This is consistent with our previous remarks. 



Acknowledgment The authors are grateful to Prof. H. D. Marta, from the 
Institute de Fisica, Facultad de Ingenieria, Montevideo, Uruguay, for 
his help in the calculation of the nuclear form factors. We thank Drs. 
Christian Beck and Monica Trotta for communicating the features of 
their new experimental results prior to their publication. This work was 
supported in part by CNPq and the MCT/FINEP/CNPq(PRONEX) 
under contract no. 41.96.0886.00. L.F.C. acknowledges partial support 
from the FAPERJ, and M.S.H. and W.H.Z.C. acknowledge support 
from the FAPESP. 

References 

[1] C.A. Bertulani, L.F. Canto and M.S. Hussein, Phys. Rep. 226, 281 
(1993); C.A. Bertulani, M.S. Hussein and G. Miinzenberg, "Physics of 
Radioactive Beams" (Nova Science, New York, 2001); M. S. Hussein, 
L.F. Canto and R. Donangelo, Nucl.Phys. A722, 321c (2003) 

[2] S. Nakayama et al. Phys. Rev. Letters 85, 262 (2002). 

[3] T. Aumann et al. , Phys.Rev.C59, 1252 (1999). 

[4] M. Trotta et a/., Phys. Rev. Lett. 84, 11 (2000). 

[5] D.G. Madland and P.G.Young, Los Alamos Report No. LA7533-mb 
(1978) (unpublished). 

[6] G. D. Alkhazov et al, Phys. Rev. Lett. 78, 2313 (1997). 

[7] V.E. Viola and T. Sikkeland, Phys. Rev. 128, 767 (1962). 

[8] Yu. N. Eldyshev, V.N. Lukyanov. and Yu. S. Pol, Sov. J. Nucl. Phys. 
16, 282 (1973). 

[9] C.H. Dasso, S. Landowne and A. Winther, Nucl. Phys. A432, 495 
(1985). 



11 



[10] K. Hagino, A. Vitturi, C.H. Dasso and S. Lenzi, Phys. Rev. C61, 037602 
(2000). 

[11] A. Diaz-Torres and I.J. Thompson, Phys. Rev. C65, 024606 (2002). 

[12] A.M.S. Breitschaft et al, Ann. of Phys. 243, 420 (1995). 

[13] G.R. Satchler, "Direct Nuclear Reactions", Oxford University Press, 
1983. 

[14] A.R. Edmonds, "Angular Momentum in Quantum Mechanics", Prince- 
ton University Press, Princeton, New Jersey, 1974. 

[15] E. Switkowski, R.M. Wieland and Aa. Winther, Phys. Rev. Lett. 33, 
840 (1974). 

[16] J.J. Kolata et al, Phys. Rev. Letters 81, 4580 (1998) 
[17] M. Trotta, private communication. 



12 



Figure captions 



Figure 1: ^He+ ^^^U fusion cross sections. The data of refs.[4] (solid squares) 
and ref. [7] (open squares) are compared with the calculations of the present 
work. The barrier energy is indicated by an arrow. For further details see 
the text. 

Figure 2: Coulomb coupling to the breakup channel for the ^Hc + ^^^U fusion 
cross section. Experimental results [4] are compared with a static calculation 
similar to that of figure 1, with just a scaling of the potential "^He (dashed 
line), and taking into account the fact that ^He is a halo nucleus (full hue). 

Figure 3: Coulomb and nuclear dipole form factors (a) and their ratio (b). 
See text for details. 

Figure 4: Total fusion data of ^He incident on ^^^U in comparison to the 
static calculation, including the ^He halo, of figure 2 (dashed line), a coupled 
channel calculation including only the Coulomb interaction (thin full line), 
and also including nuclear effects (thick full line). See text for details on 
these two last calculations. 



13 




Figure 1 



lo'i ' I ' I ' I ' I ' I ' I ' I ' I ' ^ 

: 6He + 238u 



t 

b 



10' r 



10^ r 



10^ 



15 




• Data of Trotta et al - 

I without static effect _= 

with static effect 



20 

^c.m. (MeV) 



25 



30 



Figure 2 




Figure 3 



1x10 E~i — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




15 20 25 30 



^c.m. (MeV) 



Figure 4 



